Proton-Nucleus Cross Section at High Energies 
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Cross sections for proton inelastic collision with different nuclei are described within the Glauber 
and multiple scattering approximations. A significant difference between approximate "Glauber" 
formula and exact calculations with a geometrical scaling assumption for very high-energy cross 
section is shown. Experimental values of proton-proton cross sections obtained using extensive air 
shower data are based on the relationship of proton-proton and respective proton-air absorption 
cross sections. According to obtained results values reported by the Akeno and Fly's Eye exper- 
imental groups are about fO% overestimated. The proper energy dependence of absorption cross 
section for collisions with air nuclei is of a great importance for studies of high energy cosmic rays 
using the Monte Carlo technique. 
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The rise of the proton-proton cross section (total, inelastic) as the interaction energy increases is an important feature 
of the strong interaction picture. The growth itself is established quite well both from theoretical and experimental 
point of view. However the question how fast do cross sections rise is discussed permanently. A definite answer is still 
lacking. Theoretical predictions agree well with one another and with accelerator data in the region where data exist 
{\fs ~ 20 ~ 2000 GeV) but they differ above taht level. Before the large hadron Collider LHC Q shifts the direct 
measurements limit to 10 ~ 14 TeV the only existing information comes from the cosmic ray extensive air shower 
(EAS) data. Two EAS experiments, Akeno |2| and Fly's Eye ||], gave estimations of proton-proton total cross section 
at about y/a » 10 4 GeV. 

The important difference between the collider and EAS proton-proton cross section measurements is that in fact 
the proton-air interactions are involved in the EAS development. Thus the value which is real measured is the cross 
section for the interactions with air nuclei. The value of proton-proton cross section is then obtained using a theory for 
nuclei interactions. In many recent papers concerning this subject one finds only a brief reference such as "calculations 
have been made in the standard Glauber formalism" or something very similar [^]||. 

The original Glauber paper Q was published over 40 years ago and it appears that since then some misunderstand- 
ings have arisen. Rather complicated equations for scattering cross sections can be simplified significantly applying 
some additional assumptions — of limited validity. It was pointed out in 1970 ||] that some expressions which are 
most frequently identified with the nuclear optical model should not be used at least for light nuclei, yet this still 
sometimes happens nowadays. 

In this paper we compare results of calculations with and without mentioned simplifications. We will shown that 
there is quite significant difference between them. The exact Glauber formalism will be used to evaluate the proton- 
proton cross section values from the Akeno and Fly's Eye data. 

The paper is organized as follows. In the next section a detail description of proton-proton scattering, the basis 
of further nuclei cross section calculations, is given. In section 3 the proton- nucleus cross section evaluation for some 
commonly used methods is given, and in section 4 the quantitative results are presented and discussed. 



II. PROTON-PROTON CROSS SECTIONS 

Introducing the impact parameter formalism cross sections can be described using one, in general complex, function 
X in the form 
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The phase shift x is related to the scattering amplitude by the two-dimensional Fourier transform 
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1 _ c *x(b) = _L f c -^t5-( t )d 2 t; 

27TI J 

S{t) = ^- J e lbt (l - e lx ^d 2 b. (2) 

Using the optical analogy one can interpret the 1 — e**^ function as a transmission coefficient for a given impact 
parameter. Considering two colliding object we can assume (for pure absorptive potential) 

X (b) = iu{b) = iK ab J d 2 b' p a (b)pb(b + W), (3) 

where ph is a particle's "opaqueness" (the matter density integrated along the collision axis). To some extend 
the hadronic matter density could be identified with the charge density measured precisely in the leptonic-scattering 
experiments. However, for high-energy hadron-hadron collisions the real part of the phase shift x c & n be considerable. 
Also, the simple interpretation given in Eq.(^|) ought to be modified — at least, by introducing some dependence on 
the interaction energy (s). Thus the general phase shift x is a two- variable complex function. Two possible ways 
of simplifying the situation have been proposed in the literature: the factorization hypothesis (FH) and geometrical 
scaling (GS). They can be expressed as 

X (s, b) = iw(b)f(s) (FH) 
X (s, b) = iu(b/b ( s )) (GS). (4) 

From the optical point of view, the FH means that the hadron is getting blacker as the energy increases while the 
GS means that it is getting bigger. 

The bulk of information about the hadron phase-shift function x comes from elastic scattering experiments: more 
precisely, from the measured differential elastic cross section. The value of the total cross section can be obtained 
from the imaginary part of the forward scattering amplitude using the optical theorem. Analysis of the clastic data 
above y/s ~ 20 GeV shows that neither assumption given in Eq.(^) is realized exactly (see e.g. Refs. However, 
a combination of the two can reproduce the data quite well. 

In this paper the form of x is assumed after Q in the form 

X (s,b) = (A(s)-MMM (5) 

[with oj defined by Eq.(||)]. This follows the original GS idea |J and differs from the known Martin formula where the 
ratio of real to imaginary part of the scattering amplitude depends on the momentum transfer. A full discussion and 
some recent references can be found in |^[ . In any case, due to a lack of both theoretical and experimental information 
about the phase shift, Eq.(|^) can be treated as a first approximation. In this paper the parametrization of A used is 



A(s) = — _ ( So = 500 GeV 2 . (6) 



0.077 ln(s/s ) rMM r% 

1 + 0.18 ln(s/s ) + .015 ln 2 (s/s ) 

For the uj energy dependence in Eq.(Eh the GS is assumed, thus 



u)(b, s) — uj(b) with 
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where s is the center of mass energy for which the detail shape of u> has been originally determined. The accurate 
data description can be found in Ref. jl0| with the "hadronic matter" distribution of the form 

Ph (b) = J dz^e- m » r (8) 

with the coefficients m h [and K ab in Eq.(^)] adjusted to the hadron-proton elastic scattering data at ^/s ~ 20 GeV. 
The energy dependence of the phase shift is thus introduced by using the elastic to total cross sections ratio change 
[A(s) function] and by the scaling of the interaction impact parameter which reproduces the increase of cross section 
as the interaction energy increases [(Tinei(s)]. 

The parametrization of (Jinci(s) has been made in the ln 2 (s) form 

cr inc i(s) = 32.4 - 1.21n(s) + 0.211n 2 (s) (9) 
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which gives the cross section in mb when s is given in GeV 2 . The deviation from the Pomeron-type power-law fit (see 
e.g. |nj) is negligible in the region where the fit can be compared with existing data. The advantage of using the In 2 
form is that we do not need to concern oruselves about the violation of the Froissart unitarity bound when using our 
formula at very high energies. 

Using Eqs.(|, § and |) we can calculate hadron-nucleon cross sections at any energy of interest from the point of 
view of EAS physics. 

The quality of proposed parametrizations is presented in Figs.fy and 0. 
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FIG. 1. Differential p-p elastic cross sections obtained using the proposed parametrization of \ f° r different energies compared 
with experimental data from FNAL, ISR and SPS [12]. 
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FIG. 2. The energy dependence of p-p total cross sections (a) and imaginary to real part of elastic forward amplitude (b) 
calculated using the proposed parametrization of A compared with experimental data from [13] . Fhe broken curve shows Block 
and Cahn fit [11] and the dotted corve Durand and Pi [14] calculation result. 



III. NUCLEON-NUCLEUS CROSS SECTION CALCULATIONS 



There are two different approaches to nuclei cross section calculations. Both give similar results, but foundations are 
quite different. Broadly, it can be said that the Glauber method is based on the superposition of scattering potentials, 
while the other method (which we shall call hereafter the multiple scattering method) works on probabilities. The 
multiple scattering approximation is widely presented in |l5| ] , mainly for the elastic scattering of different particles on 
nuclei, but the probabilistic formalism can be used to obtain cross sections of other specific processes. 

Of course we do not want to judge which method describes reality better. The comparison with experimental data 
is possible in principle, however, as will be shown, expected differences between these two approaches are comparable 
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with other effects which has to be known with appropriate precision (e.g. quasi-elastic scattering). 

The description given below is performed for the proton-nucleus interactions, but the extension to the nucleus- 
nucleus case is straightforward. 



A. Glauber approximation 

For the scattering of particle on the close many-particle system (nucleus), if each interaction is treated as a two- 
particle one, the overall phase shift for incoming wave is a sum of all the two-particle phase shifts. 

A 

XA (b, {d}) = ]T X j(b - dj) (10) 
j=i 

where {d} is a set of nucleon positions in the nucleus (dj is a position of the jth nucleon in the plane perpendicular 
to the interaction axis) . Eq. ([l(]) is the essence of the Ref . @] and in fact defines the Glauber approximation (at least 
in this paper). 

The scattering amplitude is thus given by 

S(t) = ± J e' tb d 2 b J |V({d})| 2 {l - e^W)}lpdj , (H) 

3=1 

where tp describes the wavefunction of the nucleus with nucleons distributed according to {d}. If one neglect position 
correlations of the nucleons and denotes by Qj each single nucleon density we have 

A 

l^({d})| 2 = n^K) with / ejMcpT = 1 . (12) 

3=1 

If all interactions can be described by the same phase-shift function \ then 

S(t) = ±J e itb d 2 b J li^Ul^l * b 11 },/-'d ; = 



(13) 



On the other hand, the scattering process can be treated as the single collision process with its own nuclear phase 
shift Xopt(&) 

S(t) = J e ltb {l - e ix °* tW } d?b . (14) 
The comparison with Eq. ([n]) gives 

A 

e^W = y^({d})| 3 e < ^=i x ' (b - di) JJd 2 ^ - (e^™) , (15) 



where the ( ) means the averaging over all possible configurations of nucleons {d}. To go further with the calculations 
of x pt a commonly used assumption has to be made. If we assume that the number of scattering centers (A) is large 
and the transparency of the nucleus as a whole remains constant then 

Xopt(&) = i fd 2 dp A (d) [l-e^ b - d >l . (16) 



where pa is the distribution of scattering center (nucleon) positions in the nucleus Qj)- 

When the individual nucleon opacity |1 — e lx ^\ is a very sharply peaked compared with pa then with the help of 
the optical theorem the simple formula can be found 
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Substituting (O) into (fy) the proton nucleus inelastic cross section is 
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(17) 



(18) 



This 



where the last equality holds in the large A limit [ Eq.(Q) cannot be used for A = 1 to compare results for a™ 1 ] 
result is often but not quite correctly called "the Glauber approximation" . As has been shown, the original Glauber 
assumption given in Eq.(UG) has to be supported by small nucleon sizes and a large value of A. 



B. Multiple scattering approach 



The fiinci given in Eq.([l]) can be interpreted in the probabilistic way by identifying the [1 — |e lx ' b -'| 2 ] term as the 
probability of inelastic scattering at impact parameter b. This can be extended to the interaction with nucleus in a 
straightforward way. If we denote this probability by P(b) and nucleons in a nucleus A are distributed according to 
Pa then the averaged probability of inelastic interaction with one of the nucleons is 



P A (h) 



<? A em p(b-d). 



The inelastic cross section with the whole nucleus is then 
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(19) 



(20) 



In the multiple scattering picture the point-nucleon approximation can be also introduced simplifying the cross 

f p cl then 

P A (b) = P -^aZ\ (21) 



section formula. If one puts -P(b) = S 2 (h) a™f then 



what leads to 
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(22) 



The above equation has a very similar form to Eq.(18) but the difference is also quite clear 



It is interesting to note that Eq.(22) is often called the "Glauber approximation" too (see e.g. [fL6[). 



C. Comparison with the exact Glauber formula 

Both approaches discussed above lead to slightly different formulae for an inelastic proton-nucleus cross section. To 
see how big the difference is we have calculated respective cross sections using exactly the same procedures, nucleus 
shapes and proton-proton cross sections. It is interesting to compare results for different nucleus masses and incoming 
proton energies. 

The energy dependence in the discussed approximation is introduced only via the proton-proton cross sections 
change. Respective formulae have been given in section 2 

The distributions of the nucleon position in the nuclei of the form given in the FRITIOF interaction model |l7]] 
have been used 

f {l + cxp^r-roA 1 / 3 ) /C]}' 1 for A > 16 
Qj ~ I [l + ^(S) 2 ]-p(-f) A<16 (23) 

with the parameters given there, except for A = 4 where we have used the so-called "parabolic Fermi distribution" 
J16[ and for lighter (A = 2,3) nuclei, where the simple Gaussian was used. The minimum allowable distance between 
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two nucleons was introduced (0.8 fm) which modify the "initial" (uncorrelated) distribution rather strongly, especially 
for very light nuclei, so the detail shape of the light nuclei uncorrelated density distribution is not necessary for our 
purposes. 

The validity of the nucleus description used is shown in FigJ^ where the proton-nucleus inelastic cross section data 
measured at low energies are compared with our calculations. The interaction energy is about the same that it was 
used as the reference energy so for the proton "hadronic matter density distribution" (v) estimation, so the energy 
dependence of all proton-proton cross section parameters are not involved here. 
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FIG. 3. Inelastic cross sections measured at the laboratory proton momentum of ~ 200 GeV/c compared with our calculation 
results as a function of target atomic mass. Data from Ref. [18]. 



The comparison between the simplified Glauber [Eq.(18)] and multiple scattering [Eq.([22|)] results is given in Fig. |J 
for proton-nucleus as well as for nucleus-nucleus interactions. We have chosen here the "air nucleus" as a target and 
five different projectiles which represent main components of primary cosmic ray mass spectrum. The difference is not 
very significant. The energy dependence is similar and the systematic shift is rather constant as maight be expected. 
It should be remembered that for light nuclei (especially for protons) Eq.rtl^) looses its physical basis. 
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FIG. 4. Cross section of collisions of different nuclei with the "air nucleus" calculated with simplified Glauber (full lines) and 
multiple scattering (broken lines) approximations for different interaction energies. 



The important point of this paper is to show how the point-nucleon approximation changes the results. Both 
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formulae given in Eqs. (J22|) and (|18|) are in fact in agreement with the factorization hypothesis for individual nucleon- 
nucleon \ function [Eq.(0)], according to which nucleons get blacker as the interaction energy increases. Our analysis, 
presented in section 2, strongly favoured the geometrical scaling which treats nucleons as getting bigger. Nucleus 
profiles obtained using exact Glauber formula [Eq.jT^)] differ from the Aa™ 1 ^ suggested by Eq.(p2|) [Aer*° t£ |- in 
Eq.(|l8|)]. The difference can be seen in Fig.|[ 




FIG. 5. Nitrogen nucleus profile functions (x p -n) obtained using the exact Glauber formula with geometrical scaling (full 
curves) and factorization hypothesis (simplified Glauber and multiple scattering) (broken curves) for different interaction energy 
(per proton-nucleon collision). 

A significant change of the nucleus size have to influence the value of inelastic cross section. Fig.|| shows the 
change in inelastic cross section of proton-nucleus with the interaction energy calculated using geometrically scaled 
and factorized nucleus profiles. 
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FIG. 6. Cross section of collisions of different nuclei with the "air nucleus" calculated using the exact Glauber formula with 
geometrical scaling (full curves) and simplified formulas (broken curves) as a function of interaction energy. 



As can be seen, the difference at very high energies is remarkable. 



IV. PROTON-PROTON CROSS SECTION FROM COSMIC RAY DATA. 



Results presented in the previous section indicate the importance of re-examination of the proton-proton cross 
section estimation based on proton-air data measured in EAS experiment. 
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The conversion from proton-air to proton-proton cross section presented in Figjj is obtained using the exact 
Glauber formalism. 
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FIG. 7. The relationship between inelastic proton-air cross section and the value of proton-proton cross section (inelastic 
- full curve and total - broken curve) calculated using exact Glauber formula with geometrical scaling. Full horizontal line 
represents the value measured in Fly's Eye [3] experiment (dashed area shows la bounds). 



The original Fly's Eye estimation of proton-proton total cross section given in is 120 mb, yet according to results 
given in Fig.0 it is 109 mb. The same procedure was applied to the Akeno data and all their proton-proton cross 
sections also appear to be about 10% overestimated. 

High-energy points plotted in Fig.^(a) which are often reproduced in the literature are taken from the original 
works. They agree quite well with the two phenomenological descriptions plotted there. The approximation of 
the cross section rise used in this paper [Eq.([|)] agrees with the Block and Cahn and with the Durand and Pi 
parametrizations in the low-energy region, but falls below the original EAS proton-proton data points. 

In Fig.|| the calculated proton-air cross section energy dependence is given. The full curve represents result obtained 
using the exact Glauber formula [Eq.([l5|)] and proton-proton phase shift x function described in section 2 [with the 
proton-proton cross section rise defined by Eq. (^)and shown in Fig.|| by the full curve] . The outcome of the "simplified 
Glauber" approach [Eq.(|l7|)] is given for a comparison by the broken curve. 
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FIG. 8. Inelastic proton-air cross sections calculated using exact Glauber formula with geometrical scaling as a function of 
interaction energy (per nucleus-nucleus interaction). Experimental points are from Akeno and Fly's Eye experiments (squares 
and the circle, respectively). Results of calculations with the "simplified Glauber" formula with the same proton-proton cross 
section energy dependence are given by the broken curve. The default CORSIKA proton-air cross section is shown for a 
comparison by the dotted curve. 
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As can be seen, the proposed cross section rise gives very good agreement with EAS measurements at about 
y/a ~ 10 4 GeV. 

Comparison of our results with outcomes from some of the most popular shower development codes (see e.g. [l9| ] 
with different models used there) shows that energy dependence of nucleus-air cross sections used there are rather 
flat, which perhaps suggests that a "simplified Glauber" have been used in cross section calculations (see the dotted 
curve in Fig.||). The relatively correct, fast, rise of the SIBYLL |2(J model proton-air cross section is obtained with 
"simplified Glauber" formalism, but with a quite extraordinary rise of the total proton-proton cross section at high 
energies (reaching about 150 mb at \/s ~ 3 • 10 4 GeV). 



V. SUMMARY 



We have shown that the geometrical scaling hypothesis with the exact Glauber formalism gives the value of proton- 
proton total cross section at about 30 TeV slightly (10%) smaller than that reported in original Fly's Eye and Akcno 
papers. The fit for <r™ el given in Eq.(||) leads to 0^°* shown in Fig.||(a) by the full curve and correctly reproduces 
Cp-air cross section measured by EAS experiments as it is shown in Fig.|[ 

The rise of the a A-air predicted by the geometrical scaling hypothesis with the exact Glauber formalism is sig- 
nificantly faster than that which can be obtained using simplified formulas (Fig.||). This can change the physical 
conclusions based on Monte Carlo simulations of the EAS development at very high energies. 
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